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I. INTRODUCTION 

The fracture functions approach [l| allows, within the framework of perturbative Quantum-Chromodynamics, 
to describe low-pj_ target fragmentation along with semi-hard radiation produced by the active parton parte- 
cipating the hard, hadron-initiated, scattering. Fracture functions are represented as non-perturbative distri- 
butions that, once properly generalized to describe multi-particles production, allow an improved description 
of QCD initial state radiation. In order to specify the connection of those distributions with the standard 
perturbative QCD approach to the dynamics of initial state, let us consider a fast incoming hadron entering 
into a high energy reaction. In a QCD-improved parton model approach, the hadron can be represented as a 
cloud of quasi-real partons which may fluctuate into far off-shell states. If no hard scattering takes place, the 
virtuality of every single space-like parton would lead the complete cascade to reassamble. If instead a hard 
scattering occurs, the virtual probe, being able to resolve parton fluctuations in the hadron up to scales of the 
order 0{1/Q^), knocks out a parton from the incoming hadron. As a result, the remaining off-shell partons in 
the initial state cascade cannot reassamble themselves any longer and therefore materialize by starting a par- 
tonic initial state shower. This radiation piles up in the final state with the one generated by the struck parton 
time-like decay in the region of 0{Q^) transverse momenta. While the latter, due to its detailed knowledge 
from e~^e~ process, is under better theoretical control, the former is at present still under study, especially in 
the kinematical regime of present and forthcoming hadron colliders. 

In a standard QCD approach to hard processes, initial state radiation is inclusively summed over by using 
DGLAP evolution equations ■ The complementary approach in which initial state radiation is unintegrated, 
or only partially integrated, gives therefore access to a more detailed inspection of space-like parton dynamics 
at high energy. The aim of this work is thus to introduce multi-particle semi-inclusive fracture functions, which 
can describe the unintegrated initial state radiation, and to give the corresponding QCD evolution equations. 

We use the general theoretical framework of Jet Calculus, originally proposed in Ref. fs*] in order to describe 
jet fragmentation. Within this framework, a highly virtual time-like parton, generated in a hard process, 
degrades its virtuality emitting a tree-like cascade of coloured quanta of lower mass until the non-perturbative 
limit is reached and hadronization takes place. Within a leading logarithmic approximation, n-particle cross- 
sections are given, in the portion of phase space specified by the jet and neglecting for simplicity correlated 
fragmentations, in terms of single-hadron fragmentation functions. The tree-like structure of the cascade is 
dictated by leading logarithmic enhancements to the cross-sections given by planar diagrams in the coUinear 
limit. Such enhancements are resummed at all order by DGLAP-like evolution equations 
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In the following we apply the same techniques in the space-like case as appropriate for the description of 
initial state radiation. In order to simplify the notation it is convenient to replace the evolution variable 
with 



5^ = TT^ln 



At LL accuracy the running of is taken into account at one loop by using 

where is the renormalization scale, Ag^^ is the QCD infrared scale, /3o = {\\Ca ~ 2ny )/127r is the one loop 
/3- function parameter and Ca and are the number of colours and flavours respectively. We then define the 
variable y 



■In 



(3) 



where Qq and stand for two arbitrary perturbative scales. We next introduce the perturbative parton-to- 
parton evolution function Ef (x, y) which expresses the probability of finding a parton j at the scale with 
a momentum fraction x of the parent parton i at the scale Qq. The function E^ (x, y) satisfies a DGLAP-type 
evolution equation (3)] 

|.^|,„,.£^^p.(„,i,,?(f.,). (4) 

where (u) are the space-like splitting functions ^2] ^nd a sum over the parton index k is understood. The 
daughter partons radiated in the evolution process of the active parton, i.e. the virtual cascade of Sec. IH are 
inclusively summed by eq. The evolution equations can be itcratively solved by using the initial condition 

El{x,y)\y^^^5l5{\^x). (5) 

The function Ef{x,y) resums large collinear logarithms of the type In"' {Q^ / Qq) . Neglecting for simplicity 
the running of as, its expansion at first order in fact reads 

Eiix, y) ^ Eiix, Ql Q') SjSil ~ x) + ^/f" (x) In ^ + O(a^). (6) 
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The E's function in LLA therefore well describe the emission of partons that are soft or close to the directions 
of the emitting parton, while large angle emissions should be included via exact matrix elements. Furthermore, 
the i?'s function satisfies the following renormalization group property: 

Eiix,Ql,Q^)^ j —E^ix/w,QlQf)El{w,QlQ'), (7) 

where the scale is such that Qq < Qi < Q^. This property is easily verified once an ag-expansion is 
performed on both sides of eq. ^ . 



II. EVOLUTION EQUATIONS FOR ORDINARY FRACTURE FUNCTIONS 

Ordinary fracture functions have been introduced in Ref. [l| in order to give a QCD-based description of 
semi-inclusive Deep Inelastic Scattering in the target fragmentation region. M^^p{x, z,Q^) represents the 

conditional probability of finding at a given scale Q^ a parton i with momentum fraction x of the incoming 
hadron momentum P while a hadron h with momentum fraction z is detected. All-order factorizazion of collinear 
and soft singularities into M^^pix, z,Q^) were demonstrated in Refs. [1, Q. In Ref. @, a fixed order 0{as) 
calculation showed explicitely that the additional collinear singularities occuring when partons are produced 
in the remnant direction can be properly renormalized only introducing fracture functions. In view of the 
generalization performed in the next Section, we briefly recall the derivation of the i p{x , z , Q'^) evolution 
equation. Let us consider the DIS one-particle inclusive cross-sections I + P I' + h + X in the target 
framentation region, t = — (P — /i)^ ^ Q^, 

f — M^p(«, z, Q2) a,ix/u, Q2) . (8) 



(2) 



(A) 



(B) 



FIG. 1: The A and B term of eq.lO. The black blob denotes the parton-to-parton evolution function E. Partons indices 
are indicated and at each triple-line vertex is associated a real AP splitting functions P{u). The diagrams are at the 
amplitude square level. The top parton line enters the hard {Q^) scattering indicated by the bright blob. 



The cross-sections is expressed as a convolution of the fracture functions with the point-hke partonic cross- 
sections (T*. Fixed order calculations shows singularities, as already said, when the emitted parton is coUinear 
to the hadron remnant. The structure of singularities is rephrased in a Jet Calculus approach by summing over 
all combinations of distributions which can give the desidered final state. Therefore we may write 

Ml^pix,z,Y) = A/^_^/^(x,z,y)+M|j,^/p(x,z,r), (9) 
dnj fx \ 

—El (-, y - 2/o) M\^^ip{w, z, yo) , (10) 



Both terms are depicted in Fig. 1. The variable j/o corresponds to an arbitrary factorization scale and Fp and 
stands respectively for the structure function of a parton i inside a proton, and the fragmentation function of 
a parton / into the hadron h. P^^^u) are the real AP splitting functions Q. The A-term stems for direct hadron 
production from a fracture function, while the _B-term, sometimes called perturbative, takes into account the 
production of /i by a timelike cascade of the daughter parton I. The scale associated to this branching, being 
not observable, has been integrated over. The integration limits in eq. ([9]) are fixed by requiring momentum 
conservation. The evolution equations are easily obtained by deriving eq.Q with respect to to Y . One has 

AM^_,/^(x,z,r) ^ J^^Pnu)M\,^^p{x/u,z,Y), (12) 

1—2 

^KM/Pi^^^^^y) = j\_'^Pi{u)^%^fp{xlu^z^Y)+ (13) 

/ * -^i^Pt iu)D^ ( ^ , y) Fh i^/u, Y) . (14) 
Jj. u xyl ~ u) \x[\—u) / 

Restoring the familiar variable Q^, we obtain the evolution equation for M^^p(x, z,Q^): 

Q'^r^K/pi^^^^^Q') = ^4?^ I\ ^Pn^)Ml/p{x/u,z,Q^) + (15) 

f{u)D'^[-^^ Q^)f},{x/u,Q') . 



QQ2 h/py ' '-^ / 27r u 



2-K u x{l-u) Va:(l-M) 
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Mj^^p{x, z, Q^) does not depend on the factorization scale and satisfies its own s-channel sum rule [i[ : 

J2 f dzzMl/pix,z,Q^) = {l-x)F'p{x,Q^). (16) 
h ■' 

At the phenomenological level, fracture functions have been shown in Ref. [3] to well reproduce at the same time 
both HERA diffractive and leading proton data, thus convalidating a common perturbative QCD approach to 
these particular classes of semi-inclusive processes. 



III. EVOLUTION EQUATIONS FOR DI-HADRONS FRACTURE FUNCTIONS 

Let us consider a double-inclusive Deep Inelastic Scattering process I + P ^ I' + hi + h2 + X where two 
detected hadrons, hi and h2, have both ^^=1.2 = — (P — /li)^ ^ Q^. In analogy with eq. ([8]) we may write the 
corresponding double-inclusive cross sections as 

-^K^M/P^''^ ^1' -^2, Q') <7j [X/U, . (17) 

M^^ fi^lp{x, zi, Z2, Q^) gives the conditional probability of finding an active quark i with momentum fraction x 
of the incoming hadron momentum P while two secondary hadrons with fractional energy zi and Z2 arc detected. 
Evolution equations for di-hadron fracture functions A/^^ i^^^p{x, zi, Z2,Q'^) can be obtained generalizing the 
derivation outlined in the previous Section. M^^ ^^yp(a;, zi, Z2, Q^) is therefore given by the incoherent sum of 
all combination of distributions which can give hi and 7i2 in the final state, see Fig. 2: 

K,M/pi^^^i^^^^Y)= ^4,/u,/../p(^'^i'^2,i^). (18) 

X=A,B,C,D 

Explicitely the four contributions read 

^i?.i,/.2/p(^'^l'^2,>') 
^BJu,/.2/p(^'^l'^2,l') 



^C,/U,/12/P(^'^l'^2,l^) 



^DJU,/!.2/p(^'^l'^2,l') 



The v4-term is the direct convolution of a di-hadron fracture functions with the function E. The B-term involves 
one-hadron fracture functions and fragmentation functions. The third and more involved C-term is of a ladder 
type and contain two y— integration because the two hadrons are emitted by daughters of the active parton 
ordered by virtualities corresponding to yi and 1/2 ■ The D-term involves di-hadron fragmentation functions 
which obeys its own inhomogeneous-type evolution equations 0, Q. The B and C terms also contain twin 
diagrams with hi and /12 interchanged, as indicated by (/ii,zi) <-> (/i27^2)- AH the integration limits in the 
convolution integrals are fixed by requiring momentum conservation. Taking the K-derivative of all terms and 
making repeated use of the evolution equation, eq. (j4]), and of its intial condition, eq. (O, we obtain 



—-E-'- 
w ' \ w 



(w ' ^ y°)^AM,h2/p('^^ 21, 22, yo) , 



Y 



Va 



x+zi 



Ei(^,Y~y)P^\u) 
jyj u(\ ~ u) \wu / 



M 



AM/pi^^ ^^'y^^'i'i yjl^l- u) ' ^) ^ {h2,Z2), 

dx2 /-l-^i/^^ dui 



(19) 



(20) 



dy2 I dyi 



X+Z1+Z2 



dw 



x+zi ^2 Jx/X2 Wi(l Ml) 



2/ 

'2 



-P^ {U2)Ef, ( ,yi-y2]EU , F - yi 

1 — U2) \WU2 / XUi"- 



U2( 
Z2 



1X2 



w{l - U2) 



• y^)^h -..V^F'pk^, Vi) + {hi,zi) ^ (h2,Z2) , (21) 

/ 1 Va;2(l - ui) / 



dy 



P. 



dw 



X+Z1+Z2 



x/w 



w2(l — u)2 u 



%^(^,r-y). 

\wu / 



(m)F> {w, y)Dl^ i^^y : v) 

\w(l — u) w(l — u) / 



(22) 
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(A) (B) (C) (D) 



FIG. 2: The four terms of eq. (|18|l . The black blob denotes the parton-to-parton evolution function E. Partons indices 
are shown and at each triple-line vertex is associated a real AP splitting functions P{u). The diagrams are at the 
amplitude square level. The top parton line enters the hard (Q^) scattering indicated by the bright blob. 



1-^1 -Z2 



^.M'i i /P = /^^" — , " ^ Pi\u)M\^ ,p{xlu,Z2,Y)D'^'{ 

QY BMM/P I , UX(1-U) ' ^ ' / ' ^' 'I Vxl-M / 

+ (/il,Zl)^(/i2,^2)+ / —P,^u)Ml,,^^j^^/p{x/u,Zi,Z2,Y), (24) 



+ {hi,zi)^{h2,Z2)+ I ^ -^Pi{u)Ml;^^^^^^/p{x/u,z^,Z2.Y), (25) 



1-21-2:2 



^^^i., , /P - Pi\u)F},{x/u,Y)D^^'^^( 

+ / — /f(i.)Ml,^,,^,^^/p(x/z^,zi,Z2,F). (26) 



Terms containing ordinary splitting functions (u) sum up to give a homogeneous term for the evolution 
equation for Mh^ji^/p- Terms MA,h.2/p and Mb.h^/p ^^^^ ^i^^ ^Q- p4)) and eq. (|25|) combine to give 
My^^jP^ eq. We can thus finally write the evolution equation for My^^ y^^jp as: 

1-^1-^2 

+^ /* t:(I^^''W''«.',W".-'^.<3V;"(^.<?") + ('•..-'0 " + (2^) 



o / 2n .2^. («)fM^/^,Q^)^?f-"^(^f^,^^,Q^). 

ZTT Jj. U X^{\—UY \X[1—U) X{1 — U) / 

As in the case of ordinary fracture functions M^/p, eq. (|16|) . also di-hadron fracture functions My^-^ ii^/p can be 
shown to be independent on the factorization scale variable j/q a-nd satisfie s-channel sum rules: 

J2f dz2Z2Ml^j^^^p{x,zi,Z2,Q^) = {l-x-zi)Ml^/p{x,zuQ^), (28) 
J dz,z, J dz2Z2Ml,^^^^/pix,zi,Z2,Q^)^il-x)Fl,ix,Q^). (29) 

hi ,h2 
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The dihadron fracture function Mfi^j^^/p would be suitable, for example, for styding the coupled strange 
baryons/mesons production in the target fragmentation region, i.e. hi = A and /12 — K. It should be noted 
however that fragmentation functions in the strange channel are poorly known at present for one-hadron frag- 
mentation and unknown in the dihadron case. 

We wish conclude this sections by working out the general case. Consider a n-inclusive Deep Inelastic 
Scattering process I + P ^ V + hi + h2 + ■ ■ ■ hn -\- X where the n detected hadrons, hi, have all ti <Si Q^, where 
ti=i,...n = ~iP ~ and fractional energy Zi of the incoming hadron P. In analogy with eq. ([8]), we may write 
the corresponding cross sections as 

/du 
— ^/[i,.../i„/p("'^i' ■ ■ ■ <5-j(x/u,Q^). (30) 

Comparing evolution equations for Mj^/p and Mf^-^ji^/p, we obtain, by induction, the 
h /pi^T ■ ■ J ^n,Q'^) evolution equation. If M^^^p denotes the ni-hadron fracture functions and 
-Dp the n2-hadron fragmentation function, such that rii + n2 = n, we get 

Q'^^(x, zi,...,zn, Q') - /' -if Hm;/^(x^, zi,.., z„, + J2 r^h, z} ■ 
+^ /^-(^)'>f'H^Mx/.,Q^W(^,..,^,Q^). (31) 

ZTT U \X{1—U)/ \X{1—U) X[\ ~ U) / 

The inhomogeneous term in the second line of eq. pip contains a permutation Vn{h, z} over hadrons indeces 
since we have assumed to measure n distinct hadrons. The inner sum takes into account all the possible 
combinations Mn^q ® Dq which give the n-hadrons final state configurations. It can be checked that the master 
formula, eq. (j3ip . reproduces correctly the evolution equation for Mfi^ i^^j^^/p when the latter is explicitely 
calculated as done for Mj^^ i^^jp in the first part of this Section. Fair to say, this is a rather accademic exercise. 
Althought the ladder-type kinematics of LL Jet Calculus allows one to write a closed-form evolution equation 
for AIn/p, the appeareance of an increasing number of unknown distributions prevents any further analysis. 



IV. JET APPROACH TO INITIAL STATE RADIATION 



The high multiplicity problem can of course be tackled by taking advantage of the some of the properties of 
QCD radiation. As is well known, hadron activity in a given hard interaction is often coUimated in a definited 
portion of momentum space, this being a signature of the dominant coUinear branching of pQCD dynamics. 
For this reason jet cross-sections are the natural and, possibly, the most effective representation of hadronic 
final state. This approach avoids the introduction of multi-hadron distributions, which actually has caused 
the abandon of eq. pip . Perturbative calculations with an arbitrary number of partons in the final state and 
experimental jet observables can be quantitatively compared only once a common jet-algorithm is chosen and 
used on both the theoretical and experimental level. Let us sketch the jet approach to initial state radiation 
described in Ref. Q, in which the k± clustering algorithm is used. The inclusive DIS structure functions F2 

can be decomposed in terms of n-particles exclusive cross-sections i^j"' as ^10] 

00 

F2(x,g2)^^^^(")(^,Q2)_ (32) 

n=l 

If one moves from the exclusive n-particles cross-sections, eq. p2p . to the exclusive n-jets cross-sections, a 
factorized structure emerges [9| 

Ft\x,Q';Ely,^t) = E /' ^1^)4"^ |J - ycut) , (33) 

where Ucut represents the jet resolution parameter and defined in terms of an arbitrary perturbative scale E^, 
with <C Ef < Q^. In eq. ((33)) . initial state coUinear divergences are absorbed into parton distributions 
functions, Fp. The jet-coefficients R^^^ are calculable in perturbation theory and, again, are jet-algorithm 
dependent. The n-jet cross-sections can be calculated by means of fixed order calculations. The main limitation 
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of such an approach is represented by the technical difficulties of adding more and more partons in the final 
state. On the contrary, as shown in Ref. @, a leading logarithmic accurate, Jet Calculus inspired, formulation 
of n-jet cross-sections is possible. In such an approach, initial state jets, in arbitrary numbers, are accounted 
for by using a generating functional method The n-jets cross-sections are then constructed by iterating the 
block-structure G, 

G^;{U, QlQ]) ^ A^iQl plr.^^^) j^(q2^ q2) a^QI QI) (34) 



r 



FIG. 3: Representation of one jet- insertion, eq. (|34|l . The crosses represent Sudakov form factors, eq. 
parton lines vertex indicates the real branching controlled by P{u). 



\35[. The three 



for each jet insertion and schematically represented in Fig. 3. A Sudakov form factor A((5|, Q^) [ll|, defined as 



'■^0 dt f , as{t) pjk, y 



dz- 



(35) 



inhibites emissions off the struck parton lines between the two scales and Q| and thus guarantees that each 
jet is separated by the nearby jets by a inhibited-emission phase space region. Hard parton emission in eq. ([34]) 
is then controlled by real splitting functions P{u) and its subsequent time- like evolution is taken into account 
via the jet function J{Q^, fc^), defined in Ref. as 



J{Q\k')^ / dzd{z,Q',k'). 



(36) 



In eq. (|36p . the distribution d{z,Q^,k'^) expresses the probability that an initial parton, with mass Q^, decays 
into a parton with a longitudinal momentum fraction z with respect to the parent parton with virtual mass 
<C Q^. If the intermediate partons mass k^ is integrated over we get the parton fragmentation function D, 



diz,Q^,k^)=D{z,Q^). 



(37) 



In eq. ([57]) . both d and D would have only partonic indices since we are dealing with partonic jets distributions. 
Such a formulation shares with the Jet Calculus approach the iterative construction. The n-jet cross-sections 
are built using an alternate allowed and prohibited parton emissions pattern along the struck parton ordered 



virtualities, Qq < 



< 



< 



< 



< 



< Q . The main difference from the Jet Calculus approach resides 
The evolution function Ef{u,Q^,Ql), in eq. can be regarded 



in the exclusive formulation of eq. 

however as the analogous, at the inclusive level, of G'^{u,Qf,Ql) in eq. (p4|) . The former inclusively resums 
all parton emissions between the corresponding scales while the latter, instead, constrains such emissions to be 
approximately coUimated in phase space. 
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V. JET-LIKE FRACTURE FUNCTIONS 

The formalism of Ref. Q indeed can describe jets originating from hard partons decays, whose emission off the 
active parton can be controlled at the perturbative level. The beam jet is excluded from such a decription. Its 
origin is mainly of soft and kinematical nature since it results from the fragmentation of the spectator partons 
of the hadron remnants plus , eventually, semi-hard radiation coming from the evolution of the active parton. 
In the /c^-algorithm the beam jet is therefore pre-clustered and not accounted for in the n-jet cross-sections. 
In view of the importance at present at forthcoming hadron colliders of describing such part of the process, 
we propose a new semi-inclusive jet-like distribution, A^^(a;, Q^, z, f), referring to it as to a jet-like fracture 
function. Ail^{x , , z , t) expresses the probability of finding a parton i with fractional momentum x of the 
incoming hadron and virtuality Q^, while a cluster of hadrons hi is detected in a portion of phase space TZ 
defined by two variables z and t. The phase space region TZ is limited by the constraint 

7^ : U = -{P - h,f < i, to<t<Q^ . (38) 

Once the clustering procudere is performed, the variable z is obtained by summing the fractional longitudinal 
momenta of all hadrons hi satisfying the phase space constraint, eq. 



z = ^z,, h^en. (39) 

i 

In analogy with standard inclusive DIS, which makes use of parton distributions functions, we may write 

l=q.,q 

As in the inclusive case, x and are fixed by the scattered lepton kinematics. By defining the n-particle 
exclusive cross-sections, which may be obtained directly from experiments, as 

t ^ 1 d2«+2 («) 

n\dxdQ^Ul=idzmdt^ ^ ' 



we may construct the distributions in eq. (|40p by implementing the phase space constraints, eq. ((38 



1 da 



where to corresponds to beam pipe acceptance loss. At the dynamical level, the active parton described by 
increases its virtuality toward the hard vertex in a series of subsequent branchings. Futhermore, at leading 
logarithmic accuracy, each emission is strongly ordered, with the softest fcj^-emissions close to the proton 
remnant. The invariant momentum transfer constraint in eq. (|38p thus perfectly matches this leading logarithmic 
picture, including soft and semi-hard radiation into Ad]^ up to the scale t. All log-enhanced emissions above 
that scale are then resummed by using the function E, eq. ([4]). Moreover, at variance with the full inclusive 
case, the parton initiating the space-like cascade can be specified by inspecting initial state radiation itself. It 
has a fractional momentum 1 — z, where z is overall fractional momentum taken away by the hadrons with 
ti < t and the highest allowed virtuality, t, according to strong ordering. When t is chosen in the perturbative 
region we may write, according to Jet Calculus, 

Ml^{x,Q\z,t)= I —M'^{w,t,z,t) El{x/w,t,Q^). (43) 

The right hand side of eq. (j43p is depicted in Fig. 4. Convolution limits are fixed by requiring momentum 
conservation. By differentiating eq. with respect to Q^, we obtain the evolution equations for M.-'^ which 
reads 

Q'^Ml,{x,Q\z,t) = ^^ j'^^P^iu)M':,ix/u,Q^z,t). (44) 

This equation describes how the virtual photon resolves the distributions Mi^ , when the virtuality of the latter 
is varied. As already stated, the evolution equation, eq. (HH), actually resums large logarithm of the type 
Of" log"(Q^/t). In real processes, t-ordering is only partially realized. Higher order corrections produce partons 
that, even if originated by a parent parton with a hard t > t scale, could give, at the end of the time-like 
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o 



i. X. Q 



i, 1-z, t 




FIG. 4: Representation of the jet-like fracture function M\ as in eq. (|43|l . The black blob represents the inclusive 
parton emissions between virtuality t and , the bright blob on top the hard scattering process. Also indicated are the 
space- like parton momentum fractions. The forward radiation off A^<, satisfying the constraint eq. (|38p . is shown. 



shower, final state hadrons with < t. The clustering procedure of course does not distinguish the origin of 
such hadrons, nor the virtuality of the parent parton emitted off the space-like chain. All of them are however 
included in A^^, according to eq. p8|l . The description of A1< therefore becomes increasingly reliable as 
much as the accuracy in space-like partonic shower is enhanced. This, of course, can be achieved by inserting 
appropriate higher loop splitting functions 

Pl{u)^Pl^'\u) + ^Pl^'\u)^... (45) 
zvr 

in the evolution equations, eq. (|44)) . We also note that the latter is formally equivalent to the one for one-particle 
inclusive extended fracture functions of Ref. p^ . This formal equivalence is expected since we may consider 
the hadrons contained in 7?. as a pseudo-particle specified by fractional longitudinal momentum z and invariant 
momentum transfer t. At the experimental level, the proposed distributions have been already adopted in 
diffractive DIS measuraments at HERA, see Refs. [13, When the diffractive event is tagged by observing 
a rapidity gap, the unmeasured low-mass proton excitation and eventually soft p± fragments are taken into 
account by using an integrated distribution as A^<. Of particular interest is the issue concerning the AAi^ 
factorization properties and we will discuss it in the remaining part of this Section. Factorization of ordinary 
fracture functions were demonstrated in DIS in the single^article case in Refs. [1, Since is actually 
more inclusive than ordinary extended fracture functions [13|, for which factorization holds, we does not expect 
any factorization breaking effect in DIS. As a result, as for standard parton distribution in inclusive processes, 
factorization guarantees that once the beam jet and the forward radiation are measured in a given experiment 
and assigned to , such a distribution can be used in a different experiment solving tVI^ evolution equations. 
Morover factorization allows a generalization of eq. to include also the beam jet in the n-jet cross-sections, 
by simply substituting a parton distribution with a jet-like fracture function: 

F^^+^\x,Q^:Ely^^^^Y. f '^^ l\ — M^^{x/u,Q^ z,y,^tEf) R^^^ (u, as, Vent) ■ (46) 

In eq. (|46|) . the fraction z of the beam jet, being not measured, is integrated over. The scale t is set to t ~ Vcut^t 
and A^!j should be properly evolved according to eq. (|^^ before being inserted in eq. pS)) . Turning now to 
hadron-hadron collisions, consider a semi-inclusive Drell-Yan type process: 

P + p-^Ci+C2+7*^ i+i- + X . (47) 

The invariant mass of the lepton-pair provides the perturbative hard scale and the leptonic final state allows a 
clean inspection of QCD initial state radiation. Ci and C2 are two hadronic forward clusters, containing each 
a beam remnant jet and associated radiation, and are defined in the phase space regions and 7^2, eq- ([38|) . 
limited by ti and t2 ■ If we assume for a while that factorization holds, the cross-sections for the process can be 
written as 



da 



DY 



i dxidx2(Ml^{xi,Q^,zi,ti)Mii{x2,Q^,Z2,t2)+i^j^6{s-xiX2Q^), (48) 



dtidt2dQ^dzidz2 

'■,3=1,1 

in complete analogy with the inclusive Drell-Yan case, for which instead factorization has been proven in 
Refs. [3, II3]- EIQ- pictorially represented in Fig. 5, does not take into account possible multiple hard 

interaction, modelled recently in Ref. |18| , nor any model for the soft remnant-remnant interaction, which can 
be instead modelled as in Ref. [l^. Eq. (|48|) can be used to measure the strenght of factorization breaking 
effects in semi-inclusive hadron-hadron collisions since factorization properties of the cross sections could be 
connected with the observed forward radiation pattern. 
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FIG. 5: Representation of a semi-inclusive Drell-Yan process, as in eq. (|48|l . The dashed hne indicates the virtual boson. 



VI. CONCLUSIONS 

In the present work we have developed an alternative QCD-based approach to initial state radiation in hard, 
hadron-initiated process by using the Jet Calculus and Fracture Function formalism. We propose novel jet-like 
Fracture Functions which depend explicitely on a minimal set of variables describing forward radiation and target 
remnants. The lack of knowledge on soft hadronic dynamics compels the use of non-perturbative distributions. 
As in the inclusive case, the Q^-evolution of the proposed semi-inclusive multi-particle distributions can be 
predicted in QCD. As a result, an explicit description of the QCD initial state radiation becomes possible. In 
this paper we limit ourselves to results within leading logarithmic approximation. We do not foresee however 
any serious limitation in implementing higher order corrections and eventually coherence effects in the formalism 
presented here. Let us add that, since the issue of QCD factorization in hard hadron-hadron reactions is closely 
connected with the pattern of soft and semi-hard initial state radiation, this novel approach may result in a 
useful framework to reconsider it. This subject is at the moment under scrutiny as well as further possible 
applications of the formalism to the dynamics of the minimum bias and of the underlying event at hadron 
colliders. 
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